Digital PM Spectra by Transform 
Techniques 
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(Manuscript received September 10, 1968) 

We derive an expression for the -power spectrum of a class oj nonstationary 
processes with periodicity in the two-dimensional autocorrelation junction 
such that R(t 1 , t 2 ) = R(ti + T, t 2 + T). Such a class includes many oj the 
common digital signals. The method of derivation is based on the double 
Fourier transform which relates the spectrum of any signal to its autocor- 
relation function. This points to a very simple method of finding digital 
spectra. 

The results are applied to derive the general expression for the power 

spectrum of a wave phase modulated with a pulse stream ^ a n g(t — nT). 

The only restrictions on the pulse stream are that the a n 's are independent 
and have identical probability distributions, and g(t) is integrable and of 
finite length. 

I. INTRODUCTION 

The power spectrum of a signal x{t) can be defined in many ways. 
Every definition, however, has to yield some measure of the expected 
power at the output of a narrow bandpass filter, as a function of the 
center frequency of the filter. 

If x{t) is deterministic then the square of the magnitude of its 
Fourier transform represents energy density as a function of fre- 
quency. If the signal has finite length the energy is finite and we can 
define the power density to be the energy density divided by the 
length of the signal. If the signal has infinite length the energy may 
be infinite, but for realizable signals we can still define the power 
spectrum by operating on a finite time interval and finding the limit 
as the interval approaches infinity. This limit may include a set of 
8-functions. 

If x(t) is a random signal the direct way of defining the power spectrum 
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is to find the Fourier transform of a sample function x (t) on a finite time 
interval T , take the magnitude square, divide by T , average over all 
possible x (t), and finally take the limit as T — > » .* 

If x(t) is stationary, the power spectrum is proportional to the 
Fourier transform of the autocorrelation function. This is a very 
useful property which often simplifies the task of finding the power 
spectrum. 

The use of transform techniques can be extended to nonstationary 
processes by means of a double Fourier transform. 2 There is a very 
simple relationship between the double Fourier transform of the auto- 
correlation function R(t lt t 2 ) of x{t) and the expected energy (or 
power) as a function of frequency. Through the proper definitions 
this technique includes both stationary and deterministic processes 
as special cases. 

One would not expect that the term power spectrum would have 
much meaning in the general case of nonstationary processes. For 
one thing it would require infinite observation time to measure the 
spectrum. For some special classes of nonstationary processes we can 
talk about power spectra. For instance if the autocorrelation func- 
tion is a function of the time difference and only slowly varying with 
time we can talk about locally stationary processes. 

Another class of signals where power spectrum has a meaning is 
where periodicity in R(t lf t 2 ) exists. We will study signals for which 

R(h , U) = R(h + T,t 2 + T). (1) 

We will show that for this class of signals the transform technique 
can be used to derive a simple expression, equation (26) , for the power 
spectrum. It is believed that this method of arriving at the power 
spectrum is simpler than the direct way used by Anderson and Salz 
in their treatment of digital FM spectra. 3 

Equation (26) is given in such a form that it can easily be used 
for many of the common digital signals. It should be especially use- 
ful when the digital pulses are overlapping. We apply it to digital 
PM which was not included in Ref . 3. 

H. GENERAL CONSIDERATIONS 

A signal x(t) has the autocorrelationf 

RfoiU) = E[x(U)-x*(h)]. (2) 



t The symbol* denotes complex conjugation. 
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We definef 

/ R{t x , t 2 ) exp (— juitt + ju 2 t 2 ) dt x dt 2 . (3) 

• oo ** — 00 

Papoulis has derived a number of properties for the transform pair 
(see Chapter 12 of Ref. 2). 
Define 



and 



P (r) = f R(t + r, t) dt (4) 

J-oo 

W(») = ± /_" p(r)e-'" T dr. (5) 



W(w) is the average energy spectrum and it can be shown that 

PP(«) «~r(«,«). (6) 

lix(t) has infinite length we can define the average power spectrum as 

S(») = lim EM (7) 

where TP„(w) is the energy spectrum of x(t) taken over an interval of 
length T . 

Combining equations (2) , (3) , (6) , and (7) we have 

/1\/2 nTa/2 
/ E[x(h)'X*(t 2 )] exp [-jcofa - t 2 )] dk dt 2 

(8) 
If we change the order of integration and expectation and observe 
that the double integral is separable and that one integral is the con- 
jugate of the other we get 



Tl r To/2 l 2 ~l 

E\ \ x(t)e-'" 1 dt\ \ 

LI J -To/2 I J 



S(«) = lim t -'"- 7 % T LJ (9) 

which is the definition of power spectrum given by Rice. 1 
A large class of digital signals has the property 

Rit, , t 2 ) = R{1, + T,t 2 + T). (10) 

Papoulis has shown that in this case T(coi , o> 2 ) consists of line masses 
t Notice the signs in the exponential. 
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in the (co, , w 2 ) -plane. The line mass on wj = a> 2 gives the average power 
spectrum S(u>) such that 

r(«! , a> 2 ) = r r ( Wl , w a ) + (2 7 r) 2 5(a Jl ) «(«, - a> 2 ) (11) 

where r r (u>i , o> 2 ) has no line masses on co, = co 2 . From the definition of 
r(co, , w 2 ), equation (3), we get 

/CO *CO 
/ /?(/, , / 2 ) exp (-jWi + iw 2 / 2 ) d/, dt 2 . (12) 
-co ■* — CO 

Let 

k - * + t (13) 

fe = / (14) 

then 



r(w, ,«»)=[ f #(* + r, 



•exp [— jojit — j(o), — w 2 )/] dr d/. (15) 
We divide the f-axis in intervals of length T and get 

r( Wl ,co 2 ) = E / / R(t+ r,t) 

•exp [— jcoiT — j(o), — u 2 )t] dr dt. (16) 

From equation (10) it follows that we have periodicity in t; in each 
interval t e {qT, (q + 1) T] we let 

and equation (16) becomes 

r(«, , Wa ) = J) | f ft(* + r, exp (-Mr) 

•exp [-;(«. - «»)(/ + gT)] rfr d< (17) 
or 

r(w, , Wa) - / / B(/ + r, exp [-Mr - i(cu, - u 2 )t] dr dt 

T £<-—>- 7 s )" (18) 

The set of 5-f unctions gives us the line masses mentioned earlier. S(o>) is 
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the mass on <*>! = o> 2 , and from equations (18) and (11) 

S(co) = w Co £ „ ^ (/ + T ' t)e ~ iaT dT dL (19) 

We return to the original variables t x and U which are absolute times. 
Using equations (13) and (14) we get 

8(a) = 7T7fi [ [ Rfo , k) exp [-**(/, - / 2 )] dt, dh • (20) 
We segment the ii-axis in intervals of length T 

1 oo M«+nr i>t 

«(«) - cTr Z / / fl& . « 

Z7Ti „--«, J 1,-aT •>t,-0 

•exp [-jwtf, - / 2 )] dt 2 dU . (21) 
The kth term of the sum is 

&(«) = Trr; f [ R(U . O exp [-jtoft - / 2 )] d/, d/ 2 . (22) 

Z7Tl Jtr-kT •'1,-0 

From the definition of autocorrelation of equation (2) : 

E(ti,fe) = R*{t 2 ,h). (23) 

Substitute equation (23) into equation (22) and change the variables 
such that 

tx -» ti - fcT and t 2 -^ t, - kT. 

Because B(£i, U) is periodic as shown by equation (10) it follows 
that 

SM = 7T7F [ [ R *& • **) ex P \Mh ~ <*)] dt* dh (24) 

ZtI J tl --kT •'1,-0 

or 

S»(«) = flf- w («). (25) 

We can then combine terms in equation (21) where g > in pairs and 
equation (21) becomes 



S(w) = ^ |/ o / o «(/, . h) exp [-jw((, - fe)] dt, dt, 

r » «c« + i>r «r 

+ 2Re I / «(/, , / 2 ) exp [-jw(*, - / 2 )] d/, d* 2 



(26) 
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m. DIGITAL PHASE MODULATION 

Represent the signal by 

x(t) = V e cos («.< + 4>(t)) (27) 

where 

<t>(t) = S a n g(t-nT). (28) 

n- — oo 

We assume that the a n 's are independent with identical probability 
distribution and g-(0 is integrable and of finite length such that g(t) = 
outside the time interval {0, pT}. Otherwise g(t) is arbitrary. 

The Appendix shows by means of transforms that the one sided 
spectrum of x(t) is 

S(o> e + «) = y r(«, «) (29) 

where a> is the difference between the actual frequency and the carrier 
frequency. r(wi , co 2 ) is defined by equation (3) with 

fife , *») = £{exp [#fe) ~ &&)]}. (30) 

If we substitute equation (28) into equation (30) and use the fact that 
the a„'s are independent we get 

R{k , k) = ft ^[exp ija n [g(t> ~ nT) - g(t 2 - nT)]}). (31) 

n- — oo 

To show that E(ix + T, t 2 + T) = R(t lt t 2 ) let t ± -* h + T and 
t 2 -» to + T and reindex such that n -» n + 1. 

fife + F, '2 + T) 

= fl ^fexp {M. + i Wi - nT) - g(t 2 - nT)]}]. (32) 

J1-— 00 

Since we assumed identical probability distributions the periodicity 
follows. 

In order to use equation (26) we now have to evaluate R(ti, t 2 ) for 

^ U ^ °° and g ^ T. 
Since t 2 c {0, T), g(t 2 — nT) will contribute in the factors when 

-(p - 1) £ n % 0. 
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Now ti « {qT, (q + 1)7} and g(ti — nT) will contribute when 

— (p — 1) + q £ n £ q. 

In order to determine for which n's only one will contribute and for 
which n's both will contribute to the factors in equation (31) let us 
look at Fig. 1. We see that if q < p then there will be factors where 
both contribute. Equation (21) becomes 



R{U ,t 2 )= II #{exp [-jag(t 2 - nT]} 
-p+i 



• II #(exp \ja[g(ti ~ nT) - g(t 2 - nT)}}) 

-p+8+l 

•ilElexplJagfo -nT)}}. (33) 

i 

In the first product we let n -» n — p and in the second n -» n — p + 
q. We also let t 1 —^t 1 + qT and this part of equation (26) becomes 

\ Re{l!e-'- r [ T I' fl [^(exp \jag\t, + (q - n)T]}) 

Tl ^, a _l Jo Jo n-1 

•tf(exp {-jag[t 2 + (p - n)T]})] 

■ f[ E[exp (ja{g[t: + (p - n)T] - g[t 2 + (p - q - n)T}})] 

n-l 

•expf-M/, - t 2 )}dhdt 2 y (34) 

The first double integral in equation (26) we get from equation (33) 
with 5 = 0. Let n — » — n, 

-^ f f n#(exp \M9(ti + nT) - g(t 2 + nT)]}) 

Z7TI Jo Jo n-0 

•exp [-jco(h - t 2 )] dt, dt 2 . (35) 

We have yet to find the terms where p ^ q ^ °o . If we look at Fig. 1, we 
see that 

R(ti . **) - IT ^{exp \-jag(t 2 - nT)}} 
-p+i 

• fl Elexpijagih - nT)]} . (36) 
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g(t 2 -nT) 



-*l 



t i i i i i i — i i i i — n 



i I I i i i i t i i 

-p -p+q 6 q. 

Fig. 1 — Contributions to R(ti, t a ). 

In the first product let n — » — n, and in the second let n — > — n + q, and 
finally i x — > ti + qT. Substituting the result into terms in equation (26) 
where q ^ p we get 



Re 






-4 /" II #(exp [jag(t + nT)]} .«-'-' dt 

= "4 I f ' fi^lexp [w(/ + nT)]}<f"" 
Re [£«-'-'- £e-'" 

L o-i o-i 



tf/ 



(37) 



which becomes 

• T p-1 



U rn^{expfja^ + nT)]}e-'" u * 

'■* I >>Q n-0 



d/ 



sinfe^ 



r A J 27rml . 1 f (p - 1)<*T ~\ °^ V 2 / 

? „£. 5 L W - "H + 2 ~ cos L 2 — J" " 



2 J 



(38) 



Substituting equations (35), (34), and (38) into equations (29) and (26) 



V 2 
2 


1 
2tT 



S(o c + «) = Y <^T / / ft tf(exp {jofoft +nT) - ff (/,+nT)]|) 



I r r p_1 
exp [-jw(A - t 2 )] dtx dt 2 +\ Il-^lexp [jag{t + nT)]}e- iat 

I •'O n=0 



dt 



1 - 2 • cos 



r -i sin 

Sin 



I. fel\ 



\ 2 / 



2/ • 



+ 2/^2>-' au: 



DIGITAL PM SPECTRA 



405 



'[ [ II [E(ew \jag[U + (g - n)T]\) 

■A*(exp \-jag[t 2 + (p - n)T]})] H tf [exp (ja{g[h + (p - n)T] 

- ^[Ai + (P - g - n)T]})] exp [-ju(t t - t 2 )) dh dt 2 J 

+ aE f il^{exp [Ml + nT)]} 

■* m--oo I •'0 n-0 



•exp i -; 



. 2irint 



T 



(II 



S\u - 



2irm 
T 



(39) 



IV. EXAMPLE OF PM SPECTRA 



A complete evaluation of equation (39) is not attempted here. We 
compute S(w e + co) for the case when g{t) is a rectangular pulse of dura- 
tion t < T, and height 1. We assume that a has r levels which are 
equidistant and equally probable such that 

A: - 1 



O w = a + 

r 

P{a ik) ] =-■ 

First we rewrite equation (39) with p = 1: 

Sfa + «) = -^ [^y/ fl r /J-^KP Ifrfoto - 9(hW 
•exp [— jw(/i — to)] dti dt 2 
+ ^1 J* E[e ia ° il) ]exp[-jwt]dt 

•(! + -¥} -£)]■ 

The expected values in equation (42) will be 
£(exp (jtyfft) - rftf])) 

- ~ 2 exp |jja + — — ajtofo) - 0(0] 



(40) 
(41) 



(42) 



(43) 
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and 



E[e"" m ] = ) £ exp [i(a + lL y 1 a)g(£)\- 



(44) 



From equations (43) and (44) it follows that 

f [ £(exp {ja[g(h) - g(t 2 )]}) exp [-jwft - /,)] dt, 



dl 



1 ' If 5 

' i=l I •'O 



exp 



■I i k ~ 1 



g{t) — jut dt 



(45) 



since we can separate the integrals and because the integral over t 2 
is the conjugate of the integral over ti . Also 



[ E[e ,a ' w ]e~"" dt 

Jo 

f k-i Jo 



exp 



■I -L k - l 



g(t) — jut dt 



(46) 



Let us set 



f exp i(« +~~a)aW -jo/Jd* = *» (47) 



Then equation (42) becomes 

^+-)-?fe[s, i i''«i , -iSf i '*>n 

r^ 1 „ /27rm\ 5 .,' 27rm 

T 7 



+ ^Jll £r A rrV 



fi w - 



(48) 



With gU) =1 for £ e {0, t} and for £ c {t, T} we get (after some 
trigonometric manipulations) 



S(u. + u) = -r- 



r 


r . (corYl 


2 


f 




r- m 


>i 


2 

r 


sm lyJ 






1 


SUJ U 




27T7 1 


COT 

2 




1 " 


— 2 

r 


_ sin (i). 





+ E 



2 

T 


(mirT\ 
Sm \T ) 


r 2 T 2 


m-KT 
T 





M 


sm 


t) 




\2/ 




/«\ 


sm 


,^T 




\2r/J 



-:- i 1 - ^ 
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sin mir\ 1 — m 



77?7r( 1 — y 



f rT\ T 



W) sin W 1 ~ m 



sin yifrj s 
mirr 



T 



rtnr 



l-i 



sin i- 



sin i — 



cos mir + «o + 



«(r - 1) 

2?- 



]'[- ~ 



27rm "1 
T J 



(49) 



The continuous part of the spectrum is independent of «„. Only the 
random deviations from a , and the pulse length determine the con- 
tinuous part of the spectrum. The term a is related to the periodicity 
of the signal and consequently the spectrum spikes will be functions 
of a . 

There are some interesting special cases of equation (49) . If {ct/r) = 
71-2tt then the continuous part of the spectrum disappears. We would 
expect this because from equation (40) it follows that during each 
pulse a phase excursion of a plus an integer number of 2n is made. 
This is the same as modulating with a periodic signal where in one 
period we have a g (t) . 

If T -> T all spikes except the carrier spike go to zero; this gives 



S(w e + co) = -5 s 



2 



r 

2* 



sin 



wT 



2 



r 





(a\ 


sin 


Is) 




M 


sm 


!^~ 


L 


\2?7 



r 



sin 



SID 



(ll 



*(«) 



(50) 



which is a familiar result. This is independent of « because in this 
case a is just a constant phase angle added to the carrier at all times. 
It is interesting to note that if a = n-2ir, the carrier spike will also 
disappear. 



v. CONCLUSION 



There is a large class of nonstationary processes which yield sig- 
nals with periodicity in the two dimensional autocorrelation function 
such that 22 Ui + T, t« + T) = R(t u U). Such a class includes many 
of the common digital signals. 
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A simple expression has been derived for the power spectrum of 
such a signal. It was done by means of the double Fourier transform 
which relates the spectrum of any signal to its autocorrelation func- 
tion. For this class of signals the method is very powerful. 

The results were applied to get the general expression for the power 

CO 

spectrum of a wave phase modulated with a pulse stream ^Z a n(j(t ~ 

— 00 

nT). The only restrictions on the pulse stream are that the a n 's are 
independent and have identical probability distributions and, g(t) 
is integrable and of finite length. As an example a rectangular pulse of 
length t < T was considered. 

The same method can be used for digital FM to arrive at the ex- 
pression given by Anderson and Salz. 3 It can also be shown that the 
expression given here for PM goes into the one for FM given in Ref. 
3 with the accrued phase per pulse equal to 0. This corresponds to a r 
= in Ref. 3, a case which was not treated there. We get this FM 
case by substituting 



Ud 



f 9(f) di> 



for g(t) in our equation (39) where co d is a frequency deviation parameter 
defined in Ref. 3. 

APPENDIX 

Carrier Translation 
An angle modulated signal is represented by 

x (t) = V e cos [a e t + 0(0]. (51) 

The autocorrelation function of x(t) is 

BXh,t 2 ) = *[*«,) •*•&)] (52) 

which becomes 

R,(ti , k) - "^ (exp \ju.(ti + * 2 )]#{exp [#(/,) + j0(/ 2 )]} 

+ exp [-jco c (k + / 2 )]£{exp [-#(*,) - #(**)]} 

+ exp \ju,(ti ~ t 2 )]E{exp Ifrfo) - j<f>(t 2 )]} 

+ exp \-ju c (h - * 2 )]#{exp [-&(/,) + j<t>{k)\Y). (53) 
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The average energy in the two first terms will for most practical cases 
of modulation go to zero. To show this we take the double Fourier 
transform of the first term 

I\(tt, ,w 2 ) = Y^jexp Iju.ih + h)]} *5 2 (E{exp [#(/,) + J4>(k)]}) 

(54) 
or 

V* 
Ti(o)i , co 2 ) = —f- (27r)" # 5(wi — co c ) 

• 8(a> 2 + co c ) * 5 2 (E[exp [#(/,) + #&)]}) (55) 

where ff 2 is a double Fourier transform operator and * means convolution 
in both f, and t 2 . The plus sign in the second 3-function comes from the 
plus sign in front of u 2 in the definition of the double Fourier transform 
of equation (3). Equation (55) shows that we can find ri(a>! , co 2 ) by 
first finding 

ff a (tf{exp [#(*,) +J4>(t 2 )]}) 

and then move it in the (o>! , o> 2 ) -plane so that the point (0, 0) falls at 
(o e , — a) e ). (See Fig. 2.) To find the average energy we set coi = w a . 
We see then, if 

ff a (0{exp mh) + *ft)H) 

does not have any significant mass density for frequencies of the order 
of w e there will be no contribution from I\(wi , w 2 ) falling on the line 
a i = oi 2 . This will be the assumption here, that is, the modulating func- 
tions do not produce sidebands as far as oj e away from the carrier. Then 
it also follows that the second term in equation (53) will not contribute. 
Now let us set 



Thus 



If 



E{exp [#(/,) - j0(/ 2 )]} = Rih , / 2 ). (56) 

E{exp [-#(/,) + #(/,)]} = R*fr , t 2 ). (57) 



r(«, , w a ) = *JLR{U , t 2 )] (58) 

then from equation (3) it follows 

^[R*(t lf t a )] = r(-o, 2 , - Wl ). (59) 
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Taking the double Fourier transform of equation (53) then yields 

V 2 
r,(w, , co 2 ) - -j- (2tt) 2 • [5(w, - co c ) • 5(co 2 - w c ) * r(», , co 2 ) 

-f 5(a., + w e )-S(u a + co c ) * r(-co 2 , -Wx)]. (60) 

If we look at Fig. 2 we see that the convolution just means sliding 
r(w! , co 2 ) and r(— co 2 , — coO along the line o> x = co 2 . We want the portion 
of r«(«i , co 2 ) that is located on the line. The sliding process takes the 
portions of V(u 1 , co 2 ) and r(— w 2 , — ",) that are already on the line 
and moves them to the points (co c , co c ) and ( — co e ,— co c ), respectively. 
Setting co, = co 2 = 03 c + co in equation (60) and performing the convolu- 
tion we get 

r,( W « + CO, C0 e + CO) = -f [r(cO, CO) + T(-CO, -CO)]. (61) 

The second part in equation (61) is just the mirror image of the first 
and the one-sided spectrum becomes 



where 



£(co c +co) = ^T(co,co) 



r(«! ,co 2 ) = ff 2 (E{exp [#(*,) - tiiU)]})- 



(62) 
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Fig. 2 — Currier translation. 
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